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ABSTRACT 

Context. Nuclear reaction rates of astrophysical applications are traditionally determined on the basis of Hauser- 
Feshbach reaction codes. These codes adopt a number of approximations that have never been tested, such as a 
simplified width fluctuation correction, the neglect of delayed or multiple-particle emission during the electromagnetic 
decay cascade, or the absence of the pre-equilibrium contribution at increasing incident energies. 

Aims. The reaction code TALYS has been recently updated to estimate the Maxwellian-averaged reaction rates that are 
of astrophysical relevance. These new developments enable the reaction rates to be calculated with increased accuracy 
and reliability and the approximations of previous codes to be investigated. 

Methods. The TALYS predictions for the thermonuclear rates of relevance to astrophysics are detailed and compared 
with those derived by widely-used codes for the same nuclear ingredients. 

Results. It is shown that TALYS predictions may differ significantly from those of previous codes, in particular for nuclei 
for which no or little nuclear data is available. The pre-equilibrium process is shown to influence the astrophysics rates 
of exotic neutron-rich nuclei significantly. For the first time, the Maxwellian-averaged (n,2n) reaction rate is calculated 
for all nuclei and its competition with the radiative capture rate is discussed. 

Conclusions. The TALYS code provides a new tool to estimate all nuclear reaction rates of relevance to astrophysics 
with improved accuracy and reliability. 

Key words. Nuclear Reactions - Nucleosynthesis 



1. Introduction 

Strong, weak, and electromagnetic interaction processes 
play an essential role in many different applications of nu- 
clear physics, such as reactor physics, waste incineration, 
production of radioisotopes for therapy and diagnostics, 
charged-particle beam therapy, material analysis, and nu- 
clear astrophysics. Although significant effort has been de- 
voted in the past decades to measuring reaction cross sec- 
tions, experimental data only covers a minute fraction of 
the entire data set required for such nuclear physics appli- 
cations. Reactions of interest often concern unstable or even 
exotic (neutron-rich, neutron-deficient, superheavy) species 
for which no experimental data exist. Given applications (in 
particular, nuclear astrophysics and accelerator-driven sys- 
tems) involve a large number (thousands) of unstable nuclei 
for which many different properties have to be determined. 
Finally, the energy range for which experimental data is 
available is restricted to the small range that can be stud- 
ied by present experimental setups. To fill the gaps, only 
theoretical predictions can be used. 

Reaction rates for medium or heavy nuclei, applied in 
astrophysical applications, have been calculated until now 
on the basis of the Hauser-Feshbach (HF) statistical model 
(Hauser & Feshbach 1952), adopting simplified schemes to 
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estimate the capture reaction cross section of a given tar- 
get nucleus, not only in its ground state but also in the 
different thermally populated states of the stellar plasma 
at a given temperature (Arnould 1972; Holmes et al. 1976; 
Woosley et al. 1978, Sargood 1982, Thielemann et al. 1987; 
Cowan et al. 1991; Rauscher & Thielemann 2001; Aikawa 
et al. 2005). Such schemes include a number of approxima- 
tions that have never been tested, such as the neglect of 
delayed particle emission during the electromagnetic decay 
cascade or the absence of the pre-equilibrium contribution 
at increasing incident energies or for exotic neutron-rich 
nuclei. 

In parallel, the nuclear physics community has devel- 
oped tools for specific applications, such as accelerator- 
driven systems, which can shed light on the many approxi- 
mations in nuclear astrophysical applications. One of these 
tools is the modern reaction code called TALYS (Koning 
et al. 2002, 2004, 2007). TALYS is a software for the sim- 
ulation of nuclear reactions, which includes many state-of- 
the-art nuclear models to cover all main reaction mecha- 
nisms encountered in light particle-induced nuclear reac- 
tions. TALYS provides a complete description of all reac- 
tion channels and observables and in particular takes into 
account all types of direct, pre-equilibrium, and compound 
mechanisms to estimate the total reaction probability as 
well as the competition between the various open chan- 
nels. The code is optimized for incident projectile ener- 
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gies, ranging from 1 keV up to 200 MeV on target nuclei 
with mass numbers between 12 and 339. It includes pho- 
ton, neutron, proton, deuteron, triton, 3 He, and a-particles 
as both projectiles and ejectiles, and single-particle as well 
as multi-particle emissions and fission. All experimental in- 
formation on nuclear masses, deformation, and low- lying 
states spectra is considered, whenever available (in par- 
ticular the RIPL2 database; Belgya et al. 2006). If not, 
various local and global input models have been incorpo- 
rated to represent the nuclear structure properties, opti- 
cal potentials, level densities, 7-ray strengths, and fission 
properties. The TALYS code was designed to calculate to- 
tal and partial cross sections, residual and isomer produc- 
tion cross sections, discrete and continuum 7-ray produc- 
tion cross sections, energy spectra, angular distributions, 
double-differential spectra, as well as recoil cross sections. 
We have updated TALYS to estimate reaction rates of par- 
ticular relevance to astrophysics. In contrast to other codes 
developed for astrophysical applications, TALYS avoids 
many of the approximations mentioned above; it there- 
fore provides a unique opportunity to test the robustness 
of these alternative codes. The principal advantages of the 
TALYS code and the new benefits that it provides in par- 
ticular to he calculation of astrophysical reaction rates is 
the subject of the present paper. 

In Sect. [21 we summarize the reaction formalism and the 
major nuclear ingredients of relevance in reaction cross sec- 
tion calculations. In Sect. [31 we describe the advantages of 
the TALYS code. This includes a comparison with the code 
MOST, which has been widely used in astrophysical appli- 
cations, in Sect. 13.11 The pre-equilibrium process is shown 
in Sect. 13.21 to have an important impact on the reaction 
rate. In Sect. 13.31 the (n,2n) reaction rate is calculated for 
the first time for all nuclei and its relevance is discussed. 
Finally, in Sect.[H conclusions are drawn. 

2. TALYS and nuclear reaction model 

Most nuclear astrophysical calculations adopt nuclear reac- 
tion rates evaluated by the HF model. It relies on the fun- 
damental assumption (Bohr hypothesis) that the capture 
process occurs by means of the intermediary production of a 
compound system that can reach a state of thermodynamic 
equilibrium. This compound system is then classically re- 
ferred to as the compound nucleus (CN). The formation 
of a CN occurs if the CN level density, at the excitation 
energy corresponding to the projectile incident energy, is 
sufficiently high. The reaction 1^ + a — > /' + a' represents 
the capture of a light particle a onto a target / in its state 
fi, where /1 = corresponds to the target ground state, 
leaving the residual nucleus /' and the particle or photon 
a'; TALYS estimates the corresponding cross section by the 
compound nucleus formula for the binary cross section in 
the full jls scheme, i.e. 



2J+ 1 
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In the above equations, E a , s, ttq, I, and j represent the 
projectile energy, spin, parity, orbital, and total angular 
momentum, respectively, and l ma x is the maximum pro- 
jectile lvalue. The same symbols but labelled by a prime 
correspond to the ejectile. The symbols T 7 ', II' 1 (/', IT) 
represent the spin and parity of the target (residual) nu- 
cleus, while J and II correspond to the spin and parity of 
the compound system. The initial system of projectile and 
target nucleus is designated by a — {a, s, E a , E£ , I 1 * ,11^} 
where E% corresponds to the excitation energy of the tar- 
get nucleus, a' — {a', s', E a i, E' x , II'} is similar for the 
final system of ejectile and residual nucleus. 5 T (a) = 1, if 
(— iy7roIP = II and otherwise (the same holds for the 
final system a'). In Eq. [U A is the relative motion wave 
length, T the transmission coefficient, W the width fluc- 
tuation correction factor and D comp the depletion factor 
given by 



D co mp = ^ 
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where a reac is the total reaction cross section, <j au ^^" is 
the total discrete direct cross section, and <r PE is the pre- 
equilibrium cross section. It is assumed by the TALYS code 
that direct and compound contributions can be added inco- 
herently. The formula for D comp is only applied to weakly 
coupled channels that deplete the flux, such as contribu- 
tions from the direct or pre-equilibrium processes. For de- 
formed nuclides, the effect of direct transitions to discrete 
levels is included directly in the coupled-channels scheme 
and the a dlsc -. dlrect term is omitted from Eq. [2) 

The HF formalism is valid only if the formation and de- 
cay of the CN are totally independent. This so-called Bohr 
hypothesis may not be fully satisfied, particularly in cases 
where a few strongly and many weakly absorbing channels 
are mixed. As an example, the HF equation is known to be 
invalid when applied to the elastic channel, since in that 
case the transmission coefficients for the entrance and exit 
channels are identical, and hence correlated. These correla- 
tions enhance the elastic channel and accordingly decrease 
the other open channels. To account for these deviations, a 
width fluctuation correction W is introduced into the HF 
formalism (see Eq. [!}. The different approximate expres- 
sions included in TALYS are described and discussed in 
Hilaire et al. (2003). When many competing channels are 
open, above a few MeV of projectile energy, the width fluc- 
tuation correction can be neglected. 

Each transmission coefficient T is estimated for all lev- 
els with experimentally known energy, spin, and parity. In 

that case, we simply have (T^ l , j ,(E a >)^ = T^ ltj ,{E a >). 

However, if the excitation energy E xl which is implicit in 
the definition of channel a' , corresponds to a state in the 
continuum, we have an effective transmission coefficient for 
an excitation energy bin of width AE X determined by the 
integral 



(T£ Vj ,(E a ,)) = 



■iAB~ 



dE y p{E y ,J,U)TJ{E a .) (3) 



over the level density p(E y , J, TV), at an energy E y , spin J, 
and parity n in the CN or residual nucleus. 

For increasing energy or nuclei for which the CN does 
not have time to reach thermodynamic equilibrium, di- 
rect or pre-equilibrium processes may become significant. In 
TALYS, the direct component is derived from the Distorted 
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Wave Born Approximation for spherical nuclei and the 
coupled-channels equations for deformed nuclei. The pre- 
equilibrium emission can occur after the first stage of the 
reaction but long before statistical equilibrium of the CN 
is reached. Although prc-cquilibrium processes can cover a 
sizable part of the reaction cross section for intermediate 
energies (typically above a few MeV in stable nuclei), they 
have always been neglected to estimate the reaction rates 
for astrophysical applications. Both classical and quantum- 
mechanical pre-equilibrium models exist and are included 
in TALYS. One of the most widely used pre-equilibrium 
models is the (one- or two-component) exciton model (see 
Koning &Duijvestijn 2004 for an extensive review), in which 
the nuclear state is characterized at any moment during the 
reaction by the total energy E tot and the total number of 
particles above and holes below the Fermi surface. Particles 
(p) and holes (h) are referred to as excitons. Furthermore, 
it is assumed that all possible ways of sharing the excita- 
tion energy between different particle-hole configurations 
at the same exciton number n = p + h have an equal a- 
priori probability. To monitor the evolution of the scatter- 
ing process, one merely traces the temporal development of 
the exciton number, which changes in time as a result of 
intranuclear two-body collisions. The basic starting point 
of the exciton model is a time-dependent master equation, 
which describes the probability of transitions to more and 
less complex particle-hole states as well as transitions to 
the continuum (emission). Upon integration over time, the 
energy-averaged emission spectrum is derived. These as- 
sumptions ensure that the exciton model is amenable to 
practical calculations. The disadvantage, however, is the 
introduction of a free parameter, namely the average ma- 
trix element of the residual two-body interaction, occurring 
in the transition rates between two exciton states. 

A thermodynamic equilibrium holds locally to a good 
approximation inside stellar interiors. Consequently, the en- 
ergies of both the targets and projectiles, as well as their 
relative energies E, obey a Maxwcll-Boltzmann distribution 
corresponding to the temperature T at that location (or a 
blackbody Planck spectrum for photons). In such condi- 
tions, the astrophysical rate is obtained by integrating the 
cross section given by Eq. ([TJ) over a Maxwell-Boltzmann 
distribution of energies at the given temperature T. In addi- 
tion, in hot astrophysical plasmas, the target nucleus exists 
in its ground as well as excited states. In a thermodynamic 
equilibrium situation, the relative populations of the vari- 
ous levels of nucleus 7 M with excitation energies E£ obey a 
Maxwell-Boltzmann distribution. The effective stellar rate 
per pair of particles in the entrance channel at a tempera- 
ture T, taking account of the contributions of various target 
excited states, is finally expressed in a classical notation as 



8 \V2 



N A 



N A {av)* aa ,{T) = (—) {kT)mGi{T) 



( 2/ '' + 1) 



l 2^ aa ,{E)E^ V {-^±p-)dE, (4) 



where k is the Boltzmann constant, m the reduced mass of 
the 1° + a system, Na the Avogadro number, and 



Gi(T) = (2^ + l)/(2/° + 1) exp(-E£/kT) 



the T-dependent normalized partition function. 



(5) 



Reverse reactions can also be estimated by making use 
of the reciprocity theorem (Holmes et al. 1976). In partic- 
ular, the stellar photodissociation rates for astrophysical 
applications have until now been derived entirely from ra- 
diative capture rates by the expression 



(2J° + l)(2j„ + l) 
(2/'° - 



y /kT 



G/(T) / mkT \3/2 
G r (T)\2ri?) X 



(0) 



where Q aj is the Q-value of the I°(a, 7)/' capture. We note 
that, in stellar conditions, the reaction rates for targets in 
thermal equilibrium are usually believed to obey reciprocity 
since the forward and reverse channels are expected to be 
symmetrical, in contrast to the situation that would be ex- 
perienced by targets in their ground state only (Holmes ct 
al. 1976). The total stellar photodissociation rate can also 
be determined directly from 



£ M (2/" + 1) Af 7|B) (r) expj-EZ/hT) 
V(2J/* + 1) expC-^/fcT) 



(7) 



where the photodissociation rate \? % of state /i with ex- 
citation energy E£ is given by 



A ( 7 ,a)( T ) = / cn^E,T)a^ a) {E)dE, 

<J 



(8) 



where c is the speed of light, a? a \(E) is the photodisin- 
tegration cross section at energy E estimated within the 
above-described reaction formalism, and n 7 is the stellar 
7-ray distribution described by the blackbody Planck spec- 
trum at the given temperature T. 

The uncertainties involved in any cross section calcula- 
tion are of two origins: 

— the first one is related to the reaction mechanism, i.e. the 
model of formation and de-excitation of the CN itself. 
Reaction mechanisms have compound, pre-equilibrium, 
and direct components. The compound formation is de- 
scribed by the HF theory. The pre-equilibrium and di- 
rect contributions are often neglected, in particular in 
astrophysical codes. The impact of this non-equilibrium 
component is studied in Sect. 13.21 More specifically, the 
predictions derived by the TALYS and MOST (Aikawa 
et al. 2005, Arnould & Goriely 2006) codes are compared 
in Sect. 13.11 by applying the same nuclear input mod- 
els. The added-values of TALYS are further discussed 
in Sect. [31 

— another type of uncertainty comes from the evalua- 
tion of the nuclear quantities required for calculating 
the transmission coefficients in Eqs. (|T]) - (J5J) , i.e. the 
ground and excited state properties (masses, deforma- 
tions, matter densities, excited state energies, spins, and 
parities, . . . ), nuclear level densities, 7-ray strength, op- 
tical model potential, and fission properties. When not 
available experimentally, this information has to be de- 
rived from nuclear models. 

Ideally, when dealing with nuclear astrophysics applica- 
tions, the various nuclear ingredients should be determined 
from global, universal, and microscopic models. The large 
number of nuclides involved in the modeling of some nucle- 
osynthesis mechanisms implies that global models should 
be used. On the other hand, a universal description of all 
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nuclear properties within a unique framework for all nuclei 
involved in a nuclear network ensures the essential coher- 
ence of predictions for all unknown data. Finally, a mi- 
croscopic description provided by a physically sound the- 
ory based on first principles ensures extrapolations away 
from experimentally known energy or mass regions that 
are likely to be more reliable than predictions derived from 
more or less parameterized approaches of various types and 
levels of sophistication. This is true even as new generations 
of such models are starting to be developed to compete 
with more phenomenological highly-parameterized models 
for the reproduction of experimental data (Goriely et al. 
2002, 2004, 2007; Hilaire & Goriely 2006). Only a few reac- 
tion model codes (for astrophysical applications) adopt the 
largest possible extent of global and coherent microscopic 
(or at least semi-microscopic) models. This is the case of 
MOST (Aikawa et al. 2005, Arnould & Goriely 2006) and 
now TALYS, but this is not the case of neither the origi- 
nal work of Holmes et al. (1976), nor the code SMOKER 
(Thielemann et al. 1987) nor NON-SMOKER (Rauscher & 
Thielemann 2001). We note that MOST, SMOKER and 
NON-SMOKER codes are similar (if not identical) in their 
description of the reaction mechanisms and corresponding 
approximations, but adopt different nuclear input models. 

It should be mentioned that TALYS includes the pos- 
sibility to calculate reaction rates with both different in- 
put parameters and a number of alternative local or global 
models for the various nuclear ingredients. This feature al- 
lows the possible study of not only parameter uncertainties 
(and therefore covariance matrices), but also model uncer- 
tainties. The readers are referred to Aikawa et al. (2005), 
Arnould& Goriely (2006) and Arnould et al. (2007) for a 
discussion on how nuclear models affect the reaction rates 
of relevance to nucleosynthesis. 

Finally, we note that TALYS, in a similar way to all 
the astrophysically orientated codes, does not describe the 
direct capture mechanism (Christy & Duck 1961; Satchler 
1983), which is known to play an important role in light, 
closed shell, or heavy exotic neutron-rich systems for which 
no resonant states are available. The contribution of the 
direct capture of heavy exotic neutron-rich nuclei has 
been taken into account only by a few restricted studies 
(Mathews et al. 1983; Goriely 1997, 1998; Rauscher et al. 
1998). This direct radiative capture is therefore not in- 
cluded in TALYS. The so-called direct process in TALYS 
corresponds to the elastic or inelastic excitation of the dis- 
crete levels at the highest outgoing energies, and is de- 
scribed using either the distorted wave Born approximation 
for spherical nuclei or the coupled-channel equations for de- 
formed nuclei. The low-energy direct radiative capture to 
a final discrete state, as described for example by the po- 
tential model (Goriely 1997) and involving the overlap of 
the initial scattering wave function in the entrance channel, 
the bound-state wave function and the electromagnetic op- 
erator, is not taken into account by the TALYS code, and 
is not discussed further here. We do mention however that 
a simple but powerful simulation of the direct-semidirect 
capture process is included in TALYS by the exciton model 
(Akkermans & Gruppelaar, 1985). 

3. The advantage of TALYS predictions 

Concerning the calculation of Maxwellian-averaged reac- 
tion rates of astrophysical interest, the TALYS code has 
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Fig. 1. Stellar 187 0s(n,7) 188 0s rate as a function of the 
temperature for 3 different prescriptions of the width fluc- 
tuation correction and without any. 

some clear advantages over previous codes developed for 
astrophysical applications (like MOST), such as 

— the inclusion of the pre-equilibrium reaction mechanism 
(totally neglected in other astrophysical nuclear codes), 

— the detailed description of the decay scheme, includ- 
ing the description of 7-delayed particle emission and 
the possible particle emission from all residual nuclei. 
In astrophysical codes, the particle emission following a 
primary 7 cascade is not explicitly taken into account, 
but corrected roughly by artificially reducing the pho- 
ton transmission coefficient in the numerator of Eq. [TJ 
as described by Holmes et al. (1976); 

— the inclusion of multi-particle emission (totally ne- 
glected in other astrophysical nuclear codes); 

— the inclusion of detailed width fluctuation corrections 
(Moldauer 1980, Verbaarschot 1985) in contrast to the 
approximation of Hofmann et al. (1980) applied in other 
astrophysical codes; 

— the inclusion of parity-dependent level densities (the 
parity equipartition is usually assumed in astrophysi- 
cal codes, one exception being the work of Mocelj et 
al. 2007, which however considers a spin-independent 
parity distribution); 

— the inclusion of a coupled channel description for de- 
formed nuclei, while the other astrophysical codes con- 
sider spherically equivalent optical potentials for de- 
formed targets, and 

— the inclusion of the fission channel for the compound as 
well as the residual nuclei. Fission is often neglected in 
astrophysical codes, and if included, insufficiently tested 
on experimental data and not consistently taken into 
account in the full decay scheme. 

For the calculation of thermonuclear reaction rates at 
astrophysically relevant temperatures, most approxima- 
tions assumed by previous codes were justified, since for 
neutron captures only incident energies of a few keV to 1 
MeV contribute, and for incident charged particles only en- 
ergies below or close to the Coulomb barrier come into play. 
Most of the above mentioned "added values" therefore do 
not affect the predicted rates by more than a few tens of a 
percent and remain insignificant with respect to the uncer- 
tainties in the nuclear ingredients, such as nuclear proper- 
ties, level densities, or 7-ray strengths. This concerns in par- 
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ticular the influence of effects such as multi-particle emis- 
sion, the width fluctuation correction, the parity-dependent 
level density or the coupled channel description for de- 
formed nuclei. Their impact on astrophysics applications 
may consequently be relatively limited. 

However, in some specific astrophysical applications, 
highly accurate predictions may be required. This is the 
case for the Re-Os nucleocosmochronology, which is known 
to be sensitive to the estimate of the 187 0s(n,7) 188 0s re- 
action rate and most particularly to the experimentally 
unknown contribution of the thermally populated 9.75keV 
first excited state of 187 0s. A significant amount of experi- 
mental data were gathered to constrain the determination 
of this reaction rate directly or indirectly (see e.g. Segawa et 
al. 2007 and references therein). In this specific example, to 
reconcile all of the direct and indirect constraints obtained 
experimentally, the state-of-the-art reaction model should 
be used, and in particular a proper description of the width 
fluctuations is required (Hilaire et al. 2003). The TALYS 
code enables the neutron capture rates to be estimated with 
an improved description of the width fluctuation correction 
based in particular on the exact, but very computer-time- 
consuming, Gaussian Overlap Ensemble (GOE) approach 
of Verbaarschot et al. (1985) or its excellent approximation 
due to Moldauer (1980). Figure Q] shows that the correction 
of Hofmann et al. (1975, 1980), traditionally adopted in 
all previous works related to the Re-Os chronology, differs 
by some 40% from the GOE predictions at the relevant low 
temperatures of 1 — 3 x 10 8 K. Such a difference significantly 
affects the determination of the stellar 187 0s(n,7) 188 0s re- 
action rate, and, consequently, the determination of the age 
of the Galaxy. 

In addition to this improved modeling of the width fluc- 
tuation correction, TALYS also enables us to test for the 
first time some of the approximations made in astrophysi- 
cal codes, as shown below. It is impossible to describe all 
nuclear astrophysical applications and to fully assess the 
impact the newly-derived TALYS rates can have. For this 
reason, we restrict ourselves to a comparison of the TALYS 
predictions with those derived by more traditional Hauser- 
Feshbach codes such as MOST; this demonstrates the ex- 
tent to which some of the above-listed added values can 
affect globally the reaction rate predictions and potentially 
some nucleosynthesis processes. 

We note that we do not study the fission channel. 
TALYS was extensively used in data evaluation for nuclear 
power applications and the fission channel was tested on 
experimental data in far greater detail than for other as- 
trophysical codes (see for example, Duijvestijn & Koning 
2005). In particular, it should be mentioned that TALYS 
considers deformed optical model potentials and can ac- 
count for the calculation of the penetration through a nu- 
merical non-parabolic barrier obtained from HFB calcula- 
tions using the corresponding HFB plus combinatorial NLD 
at the saddle points (Sin et al. 2007; Goriely et al. 2008). 
Finally, TALYS takes into account the possible effect of dis- 
crete transition and class II states that may affect the re- 
action cross section at low energies. However, with TALYS 
like any other code, it remains difficult to predict (with 
a satisfactory accuracy) fission cross sections on the basis 
of global input models. For this reason, a detailed study 
dedicated to fission will be presented in a future specific 
publication (Goriely et al. 2008). 
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A 

Fig. 2. Comparison of MOST (crosses) or TALYS (circles) 
radiative neutron capture rates with experimental data 
above Ca (Bao et al. 2000) at the temperature T = 3 10 8 K. 



3.1. Comparison between MOST and TALYS predictions 

The present section aims to compare TALYS predictions 
with those derived by the MOST code which has been 
widely used in astrophysical applications (Aikawa et al. 
2005; Arnould & Goriely 2006). The same nuclear in- 
put models are adopted for the comparison, namely the 
Hartree-Fock-Bogolyubov (HFB) nuclear masses (Goriely 
et al. 2007), the HFB plus combinatorial nuclear level den- 
sity (Hilaire & Goriely 2006), the HFBCS plus QRPA 7- 
ray strength (Goriely & Khan 2002), the global nucleon 
optical model potential of Koning & Delaroche (2003), and 
the a-potential of McFadden & Satchler (1966). Whenever 
available, experimental data constraining such models are 
included, although both reaction codes do not make use of 
exactly the same renormalization procedures. It should be 
noted that an accurate comparison of both reaction codes 
remains impossible due to the different recipes used either 
in the renormalization procedure or the way in which ex- 
perimental data are used explicitly in the code or implicitly 
to constrain theoretical inputs. In particular, the number of 
discrete states, as well as their spin and parity assignment 
could differ. TALYS code also renormalizes the 7 strength 
function on radiative strength at the neutron binding en- 
ergy when available. This step is not included by MOST. 

With such global inputs, both codes reproduce the ex- 
perimental neutron, proton or photo-rates with the same 
accuracy, i.e. to within roughly a factor of 2-3, as shown 
in Fig. [5] for the radiative neutron capture rates of particu- 
lar interest in s-process nucleosynthesis studies. Differences 
can however be noted. These arise partly because of the dif- 
ferent databases of experimental structure properties that 
were adopted and more particularly the way in which the 
deformation effects were treated. In MOST (but not in 
TALYS), the Fe-group nuclei are described as rotationally 
deformed nuclei rather than vibrational nuclei. This ex- 
plains the larger neutron capture rate of low-mass nuclei. 
TALYS also estimates the width fluctuation correction on 
the basis of the Moldauer (1980) prescription, which tends 
to increase (see Fig. [T]) the reaction cross section relative 
to the Hofmann et al. (1980) correction adopted in MOST. 
This effect is stronger in spherical nuclei and explains the 
tendency for TALYS to predict larger neutron rates than 
MOST for heavy-mass nuclei. 
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Although some discrepancies can be observed, both 
codes predict globally very similar radiative neutron cap- 
ture rates for the stable nuclei. This may not be the case 
for exotic neutron-rich nuclei, such as those involved in the 
r-process nucleosynthesis (for a review on the r-process, see 
Arnould et al. 2007). We show in Figs. [SHU the discrep- 
ancies that could exist between MOST and TALYS esti- 
mates of the radiative neutron capture rates for the Sn 
and Pb isotopes (Fig. [3]) and more generally for all nuclei 
with 16 < Z < 83 between the neutron and proton drip 
lines (Fig. [J). This comparison is made for a temperature 
of Tg = 1 (where Tg is the temperature expressed in 10 9 K), 
which is typical of conditions expected during the neutron 
captures by exotic neutron-rich nuclei in the r-process nu- 
cleosynthesis. As long as we deal with stable nuclei, both 
codes predict almost identical rates. However, deviations 
by a factor as large as 10 for Sn and Pb are seen, as soon 
as nuclei away from the valley of /^-stability are considered. 

These deviations can be understood in terms of 

— the sensitivity of the radiative neutron capture, espe- 
cially for nuclei with low neutron separation energy 
(S n ), to the nuclear structure properties and nuclear 
ingredients entering the reaction model. In particular, 
such a comparison emphasizes the difficulty in predict- 
ing neutron-capture cross sections when no or little ex- 
perimental data are available to constrain the nuclear 
inputs or to enter directly into the calculation. For ra- 
diative neutron capture, this is particularly the case 
for spherical targets (the HFB mass model used here, 
namely HFB-14, predicts the 142 <> A <> 163 Sn and 
223 A ^ 248 Pb n-rich isotopes to be deformed). For 
example, the 211 Pb neutron capture cross section and 
rate differ by about a factor of 7 (Fig. [3]). This devi- 
ation is essentially due to the number of experimental 
discrete levels of the compound nucleus included in each 
calculation. Different prescriptions are used to estimate 
the number of levels up to which the level scheme is be- 
lieved to be complete and, in this specific case, this has 
a direct impact on the low-energy cross section. 

— the reaction mechanism and its particular implemen- 
tation. In particular, the contribution of the pre- 
equilibrium process, totally neglected in other astro- 
physical codes, in fact significantly modifies the cross 
section of the keV-MeV region of astrophysical inter- 
est (as described below). This pre-equilibrium effect is 
also responsible for the odd-even staggering observed 
for exotic neutron-rich nuclei in Fig. [3] (as discussed in 
Sect.EE2). 



3.2. The pre-equilibrium contribution 

In the empirical picture of a nuclear reaction, it is imag- 
ined that the incident particle step-by-step creates more 
complex states of excited particles and holes in the com- 
pound system of target plus projectile and thereby grad- 
ually loses its memory of the initial energy and direction. 
Pre-equilibrium emission occurs after the first stage of the 
reaction but long before statistical equilibrium of the CN is 
achieved, leading to the well-known high-energy tails in the 
emission spectra and the smooth forward peaked angular 
distributions. The pre-equilibrium process is also respon- 
sible for the observed increase in the radiative neutron- 
capture cross section for stable nuclei at incident ener- 
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Fig. 3. Ratio of MOST to TALYS radiative neutron cap- 
ture rates for all Sn (crosses) and Pb (circles) isotopes at 
T = 10 9 K as a function of the neutron number N. 
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Fig. 4. Ratio of MOST to TALYS radiative neutron cap- 
ture rates for all nuclei with 16 < Z < 83 at T = 10 9 K as 
a function of the neutron number N. 



gies typically around 10 MeV (see e.g. upper panel of 
Fig. O, which hardly affects the Maxwellian-averaged re- 
action rates of astrophysical relevance. However, if we con- 
sider exotic neutron-rich nuclei with low neutron separa- 
tion energies, the low- level density makes it difficult for 
the nucleus to achieve complete equilibrium. In this case, 
the pre-equilibrium process starts to affect the neutron 
channel already at a few hundred keV, and consequently 
may change the astrophysical rate. An example is shown in 
Fig. [5l where it can be seen that already at 100 keV, the 
pre-equilibrium process affects the capture cross section of 
140 Sn, while for stable nuclei, such as 116 Sn, it starts to con- 
tribute only above some 7 MeV. Interestingly, for neutron- 
rich nuclei the calculated capture cross section for the keV- 
MeV range is smaller when the pre-equilibrium mechanism 
is included rather than excluded. This is because already 
at these low energies the neutron inelastic channel basically 
exhausts the total non-elastic cross section. The level den- 
sity and optical models, steering the reaction processes at 
these energies, determine that the relatively small contribu- 
tion of capture, in contrast to inelastic scattering, is lower 
for pre-equilibrium processes than for compound processes. 
We note that the pre-equilibrium model adopted here cor- 
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responds to the two-component exciton model (Koning & 
Duijvestijn 2004), which remains relatively phenomenolog- 
ical as discussed in Sect. [2j Its reliability and quantita- 
tive contribution, when applied to neutron-rich nuclei, can 
therefore be questioned. However, such an effect is found to 
be qualitatively non-negligible for exotic neutron-rich nuclei 
of small neutron separation energies, regardless of the nu- 
clear input adopted. The decrease in the reaction cross sec- 
tion for the 100 keV region affects directly the Maxwellian- 
averaged reaction rates of relevance to astrophysics as seen 
in Fig. \§\ (upper panel). At increasing temperatures, the 
Gamow energy increases and overlaps the energy region in 
which the pre-equilibrium process contributes to the reac- 
tion mechanisms. The pre-equilibrium correction is found 
to be less important for odd-N target nuclei than for even- 
N. This is essentially due to the increase in the HF cross 
section for the 0.5-2 MeV region found for even-N target 
nuclei, which results from the oddness of the compound sys- 
tem that has a larger nuclear level density and the evenness 
of the target reducing the inelastic channel competition. 

The pre-equilibrium contribution is found to decrease 
not only the radiative capture, but also the (n,2n) reac- 
tion rates of neutron-rich nuclei at temperatures around 
T ~ 2 - 4 x 10 9 K by about 50%. It affects not only the 
nucleon or a-particle capture rates but also the photoneu- 
tron rates of neutron-rich nuclei, as illustrated in Fig. [5] 
(lower panel). Since the pre-equilibrium component con- 
tributes increasingly more as energies increase, the higher 
the temperature the larger the impact on the astrophysics 
rates, as seen in Fig. [SJ The impact is also observed to be 
stronger for odd-N compound nuclei (i.e. for the neutron 
capture by an even-N target or the photoneutron emission 
of an odd-N target) . 

3.3. The multi-particle emission 

All reaction codes dedicated to astrophysical applications 
have made the simplified assumption that the contribu- 
tion of the emission of more than one particle in neutron 
or charged-particle induced reactions is negligible. In most 
astrophysical scenarios, the temperatures are indeed suffi- 
ciently low for such channels not to compete. However, it 
remained impossible to test such an approximation, because 
until now, no reaction code for astrophysical applications 
has included the multi-particle emission in their frame- 
work. By applying the TALYS code, it is possible to have a 
more robust understanding of the competing channels that 
until now have been neglected. An important one is the 
(n,2n) channel that may dominate, at relatively high tem- 
peratures, the (n,7) channel, especially for exotic neutron- 
rich nuclei. For energies above the two-neutron separation 
energies, the (n,2n) channel is open and begins to domi- 
nate very quickly, so that at a given temperature the ratio 
r n = (cv) / n ,2n) / ( av ) (n,^) becomes larger than 1. This tem- 
perature is critical in applications such as the r-process nu- 
cleosynthesis, since above such a temperature, the neutron 
capture will not lead to the production of heavier neutron- 
rich nuclei, but in contrast will start to reverse the nuclear 
flow and prevent it from reaching the neutron-rich region. 

In Fig. [Jj we show the evolution of the critical tempera- 
ture in the (N,Z)-plane. It can be seen that at temperatures 
above 3 x 10 9 K, most of the neutron-rich region cannot 
be reached by neutron captures regardless of the neutron 
density in the astrophysical environment. At Tg = 3 for 
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Fig. 5. Radiative neutron capture cross section for 116 Sn 
(top) and 140 Sn (bottom) with and without including the 
pre-equilibrium processes as a function of energy 

example, within the Cd or Sn isotopic chain, the r-process 
flow could not reach isotopes heavier than 131 Cd or 133 Sn, 
respectively. At such high temperatures, the (n,2n) channel 
represents a nuclear barrier that the r-process flow would 
not be able to cross, independently of the amplitude of the 
neutron density. The nuclear flow into the neutron-rich re- 
gion is usually believed to be stopped by photodissociation. 
However, this (7,«) barrier at a given temperature can 
always be counterbalanced by considering higher neutron 
densities, while the (n, In) competition represents a barrier 
to the (n, 7) flow independent of the neutron density. This 
nuclear barrier could significantly modify the prediction of 
the resulting r-abundance distribution, should the r-process 
develop in such a high-temperature environment. The tem- 
perature at which the r-process occurs remains however un- 
known. 

4. Conclusion 

The reaction code TALYS has undergone several develop- 
ments. These include the modification of the code to es- 
timate the Maxwellian-averaged reaction rates of partic- 
ular relevance to astrophysics, and the inclusion of mi- 
croscopic nuclear models for ground state properties, nu- 
clear level densities, 7-ray strength functions, and nucleon- 
nucleus optical potentials. These developments enable the 
reaction rates to be estimated with improved accuracy and 
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Fig. 6. Ratio of the radiative neutron capture (top) and 
photoneutron (bottom) rates obtained with and without 
including the pre-equilibrium processes for all Sn isotopes 
at 3 different temperatures Tg. 



reliability, and the approximations, assumed until now by 
other codes, to be tested. The TALYS predictions for the 
thermonuclear rates are found to be relatively similar to 
the more traditional astrophysics codes for stable nuclei, 
but major differences are obtained when considering ex- 
otic neutron-rich nuclei. In particular, the pre-equilibrium 
process, neglected in previous astrophysically-orientated re- 
action codes, is found to play an important role at low 
energies, at least for neutron-rich nuclei. TALYS can also 
estimate the multi-particle emission cross sections in pho- 
ton and particle-induced reactions that have been until ne- 
glected. More specifically, the (n,2n) channel was found to 
become important at increasing temperatures and even to 
dominate, at relatively high temperatures, the (11,7) chan- 
nel, for exotic neutron-rich nuclei. The presence of such a 
(n,2n) barrier demonstrates that the neutron captures char- 
acterizing the r-process nucleosynthesis could not occur at 
temperatures above typically T — 3 x 10 9 K. 

TALYS presents additional advantages such as the de- 
tailed description of the decay scheme, the inclusion of ac- 
curate width fluctuation corrections, the inclusion of parity- 
dependent level densities or a proper description of the de- 
formation effects. The impact of such features on reaction 
rates of relevance to astrophysics and on astrophysical ob- 
servables remains to be studied in the future. 



All of these developments can now be applied to im- 
prove the determination of reaction rates which will hope- 
fully help to place nuclear astrophysics applications on a 
safer footing. The TALYS predictions for the thermonu- 
clear rates are made available to the scientific commu- 



nity at the http://www-astro.ulb.ac.be website for neu- 
tron, proton, and a-induced reactions at temperatures 
ranging between 10 5 and 10 10 K for all nuclei between 
Z = 8 and Z = 100 lying between the proton and neu- 
tron drip lines. Th e TALYS code is publicly available at 
\http://www. talys.'eu The neutron, proton and a-particle 
capture rates for all 8 < Z < 105 nuclei lying between the 
proton and neutron drip lines are available in electronic 
form at the CDS via a nonymous ftp to cdsarc.u-strasbg.fr 
or via htt p: / / cdsweb.u-strasbg.fr/cgi-bin/qcat? J/ A+A/| 
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